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Abstract 

We show an 0*((£ + l)")-time algorithm for the channel assign- 
ment problem, where t is the maximum edge weight. This improves 
on the previous 0*((£ + 2)™)-time algorithm by Krai [4J, as well as 
algorithms for important special cases, like L(2, relabelling. For the 
latter problem, our algorithm works in 0*(3") time. The progress is 
achieved by applying the fast zeta transform in combination with the 
inclusion-exclusion principle. 



1 Introduction 

In the channel assignment problem, we are given a symmetric weight function 
w : V 2 — > N (we assume that € N). The elements of V will be called 
vertices (as w induces a graph on the vertex set V with edges corresponding 
to positive values of w). We say that w is ^-bounded when for every x,y £V 
we have w(x, y) < I. An assignment c : V — > {1, . . . , s} is called proper when 
for each pair of vertices x,y we have \c(x) — c(y)\ > w(x,y). The number s 
is called the span of c. The goal is to find a proper assignment of minimum 
span. Note that the special case when w is 1-bounded corresponds to the 
classical graph coloring problem. 

In this paper we deal with exact algorithms for the channel assign- 
ment problem. As a generalization of graph coloring, the decision version 
of channel assignment is NP-complete. It follows that the existence of a 
polynomial-time algorithm is unlikely. As a consequence, researchers began 
to study exponential-time algorithms for the channel assignment problem. 
The asymptotic efficiency of these algorithms is measured in terms of n = \V\ 
and I, we assume that i is a constant. The first non-trivial algorithm was 
proposed by McDiarmid [5j and had running time of 0(n 2 (2l + l) n ). It was 
then improved by Krai [1] to 0{n(l + 2) n ). 

Here we improve the running time further to 0*((£ + l) n fl 

We also 

show that the number of all proper assignments can be found in the same 
time bound. Note that for t = 1 the running time of our algorithm matches 



1 By O*0 we suppress polynomially bounded terms. 
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the time complexity of the currently fastest algorithm for graph coloring by 
Bjorklund, Husfeldt and Koivisto [I]. 

Our improvement is achieved by applying the fast zeta transform in com- 
bination with the inclusion-exclusion principle. The same ingredients were 
used also in a set partition problem in [I], however in our algorithm the 
fast zeta transform plays a different role. In particular, although channel 
assignment resembles a kind of set partition it does not seem to be possible 
to solve it by a direct application of the algorithm from [TJ. 

Some special cases of the channel assignment problem received particular 
attention. An important example is the L(p, g)-labeling of graphs, where 
given an undirected graph G = (V, E) one has to find an assignment c : 
V — > N such that if vertices u and v are adjacent then \c(u) — c(v)\ > p 
and if vertices u and v are at distance 2 then \c(u) — c(v)\ > q. The goal 
is to minimize max^ 6 yc(t;). Clearly, the algorithmic problem of finding 
an L(p, g)-labeling reduces in polynomial time to the max{p, g}-bounded 
channel assignment and we get an 0*((max{p, q} + l) n )-time algorithm as 
an immediate corollary from our result. In particular, it gives an 0*(3 n )- 
time algorithm for the most researched subcase of L(2, l)-labeling. This 
improves over the algorithms by Havet et al. [2] running in time 0(3.873") 
and a recent improvement of Junosza-Szaniawski and Rzazewski [3 J running 
in 0(3.562 n ) time. 

2 Deciding 

In this section we consider the decision version of the problem, i.e. for a given 
^-bounded weight function w and an integer s € N we check whether there 
is a proper assignment of span at most s. Since the case I = 1 can be solved 
in 0*((£ + l) n ) = 0*(2 n ) time as described in Q, here we assume t > 2. 

An assignment c : V — > N of span s can be seen as a tuple (Ji, I2, . . . , I s ), 
where Ij = c _1 (j) for every j = 1, . . . , s. We will relax the notion of assign- 
ment in that we will work with tuples of vertex sets (Ii, I2, ■ ■ ■ , Ik), where the 
ii's are not necessarily disjoint. We say that a tuple (Ji, I2, ■ ■ ■ , Ik) is proper, 
when for every i,j E {1, . . . , k} if x £ Ii and y £ Ij then \i — j\ > w(x, y). 

In what follows, U denotes the set of all proper tuples (Ii, . . . ,I S ) such 
that for each j = 1, . . . , s — £ + 1, the sets Ij,Ij+i, . . . , J,-^_i are pairwise 
disjoint. A tuple with the r last elements being empty sets is denoted as 
(Jl, . . . , I s - r , r )- For a subset X C V, we say that a tuple (Ji, . . . , Ij) lies 
in X when for every i = 1, . . . ,j, we have Ii C X. 

For v eV, define U v = {(h, . . . , I s ) G U : v € \Jj=i Observe, that 

Proposition 1. | O U v \ > iff there is a proper assignment of span s. 
By the inclusion-exclusion principle, if we denote U v = U — U v and 
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f\e0 U v = U, then 

«ey rev uey 

Our algorithm computes | f\<=y usm S the above formula. The rest of 
the section is devoted to computing | f\ey U v \ ^ or a gi ven se t Y . If we 
denote X = V — Y, then f\gy ^ ^ s J us t the set of tuples of U that lie in X: 

f]u~ v = {(h,...,i s )eu : h,...,I s CX}. (2) 

Our plan now is to compute the value of | f\eY U v \ using dynamic pro- 
gramming accelerated by the fast zeta transform. More precisely, for every 
i = £—1, . . . ,s and for every sequence J±, . . . , Ji-i of pairwise disjoint subsets 
of X our algorithm computes the value of 

i-(£-l) 

Ti~{Ji, ■ ■ ■ , Jt-\) = ■ ■ ■ > Ii-(£-i), Ji, . . . , Ji-i, $ s ~ i ) G U : |J IjCX}\, 

J'=i 

(3) 

that is, the number of tuples in U that lie in X and end with Ji, . . . , Ji-\ 
followed by s — i empty sets. Then, clearly, 

|f|^l = E TfiJt,...,!^). (4) 
veY j 1 ,...,j 1 ^ 1 c = x 
i^j^JinJj=tl) 

For every sequence of pairwise disjoint sets Ji,..., J^_j C X, we can 
initialize the value of r»_ 1 (Jx J • • • , Ji-i) hi polynomial time as followal 

Tfi^Jt, Je-i) = [(Ji, . . . , Je-i) is proper]. (5) 

Then the algorithm finds the values of for subsequent j = s. 
This is realized using the following formula: 

T*(Ji,..., J £ _i) = [(J 1 ,...,J e _ 1 ) is proper]- ^ J , Ji, . . . , Je-2), 

JoCXnproper( Ji,..., Je-i) 

(6) 

where proper(Ji, . . . , Ji-i) is the set of all vertices v € V\U- =1 Jj such that 
for each j = 1, ...,£— 1 and 2 E J,- we have j > w(v, x). 

Using the formula ([6]) explicitly, one can compute all the values of 
from the values of T*_ x in 0*((7 + l) |X| ) time, since there are tuples 
(Jo, ... , Ji—i) of disjoint subsets of X. Now we describe how to speed it up 
toO*(S x \). 

2 [a] is the Iverson's notation, i.e. [a] = 1 when a holds and [a] = otherwise. 
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Let S be a set and let / : 2 — > Z be a function on the lattice of all 
subsets of S. The zeta transform is an operator which transforms / to 
another function (£/) : 2 s — > Z and it is defined as follows: 

(Cf)(Q) = £ /(*)■ 

-RCQ 

A nice feature of the zeta transform is that given / (i.e. when the value of 
f(R) can be accessed in 0(1) time for any R) there is an algorithm (called 
fast zeta transform or Yates' algorithm, see [U [7]) which computes Qf (i.e. 
the values of ((f)(Q) for all subsets Q C S) using only 0(2l s l) arithmetic 
operations (additions). 

Let us come back to our algorithm. In the faster version, for each 
i = £,..., s, we iterate over all sequences of disjoint subsets J%, . . . , Jg-2 C 
X. Then the values of T^(J Xl . . . , J t _ x ) for all the 2 |X| ~ E ^i 1 Jjl sets J t _ x 

lyl W-2 17 1 

that are disjoint with J x , . . . , J^_2 are computed in O* (2 1 ^->j=i 1 jl ) time 
(that is in polynomial time per set!). To this end, we use the function 

f : 2 X ^=i Jj -> Z, where 

f(S)=T£ 1 (S,J 1 ,...,J e .. 2 ). 

We compute the function ((f) with the fast zeta transform using 0(2' 1 Zj j= 1 1 jl ) 
additions. Now, observe that by ([6]), for each Jg-i C X disjoint with 
Ji, . . . , J^_2, 

^ X (Ji, . . . , Je-i) = [(Ji, . . . , j£_ x ) is proper] • (C/) (Xnproper( Ji, . . . , Jt-x)). 

It follows that for each i = £,..., s the algorithm runs in time needed to 
perform the following number of additions: 

0( Yl 2 |x| -^=i | ^ 1 ) = 0( l) = 0(£ m ). (7) 

Jl,...,J^_2CX Ji,..., Ji-iCX 

By ([1]) it follows that the whole decision algorithm runs in time needed 
to perform 0(n(£ + l) n ) additions. The numbers being added are bounded 
by | HueY ^1 — — ^ n ' wnere the last inequality follows from the fact 
that the minimum span is upper bounded by (n— 1)£ + 1 (see e.g. [5]). Hence 
a single addition is performed in 0(n 2 £) time. 

Corollary 2. There is an algorithm which verifies whether the minimum 
span of an £-bounded instance of the channel assignment problem is bounded 
by s which uses 0*((£ + l) n ) time and 0*(£ n ) space. 
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3 Counting 



In this section we briefly describe how to modify the decision algorithm from 
Section [2] in order to make it count the number of proper assignments of 
span at most s. We follow the approach of Bjorklund et al. PQ. The trick is 
to modify the definition of U. Namely, now every tuple (ii, . . . , I s ) from U v 
additionally satisfies the following condition: 

s 

£|i>| = n. (8) 

i=i 

Observe, that then | P^eV ^| ec l ua l s the number of proper assignments of 
span at most s. Now, we add another dimension to the arrays T*: 

T^ k (J 1 ,...,J e _ 1 ) = \{(h,...J i _ {e _ 1) ,J 1 ,...,J^ 1 ,r- i )eU : J IjCX 

i=i 

i-(t-i) e-i 
and Yl \ I i\+J2\ J i\ =k }\- 

The dynamic programming algorithm from Section [2] can be easily modified 
to compute the values of TK^(Ji, • • • , Jt-i) for all i = £—1, . . . , s, k = 0, . . . , n 
and all sequences of I — 1 pairwise disjoint subsets of X. The details are left 
to the reader. 

Corollary 3. For any £ -bounded instance of the channel assignment problem 
the number of the proper assignments of span at most s can be computed in 
0*((£ + l) n ) time and 0*(£ n ) space. 

4 Finding 

In order to find the assignment itself we can solve the extended version of 
the channel assignment problem, where we are additionally given a set of 
vertices Z C V together with a function d : Z — > {1, . . . , s}. Then we are to 
check whether there exists a proper assignment c : V — > {1, . . . , s} satisfying 
c\z = d . It is not hard to modify the presented algorithm to solve the 
extended version of the problem in 0*((£+ l) n ~l-^l) time. The details are 
left to the reader. 

Now using the extended version of the channel assignment problem we 
can take any v € V \ Z and try each of the s < (n — 1)£ + 1 possible 
values of c{v ) one by one, each time using the algorithm for the extended 
channel assignment problem as a black box. When the value for v is fixed 
in a similar manner we assign the value for the other vertices of V \ Z. 
Since XylLiC^ + < (£ + l) n , the algorithm for finding an assignment 

has a multiplicative overhead of 0(n£) over the running time of the decision 
version. 
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5 Open problems 



In [6] Traxler has shown that for any constant c, the Constraint Satisfac- 
tion Problem (CSP) has no 0(c ra )-time algorithm, assuming the Exponential 
Time Hypothesis (ETH). More precisely, he shows that ETH implies that 
CSP requires of 1 ^ time, where d is the domain size. On the other hand, 
graph coloring, which is a variant of CSP with unbounded domain, admits a 
0*(2 n )-time algorithm. The channel assignment problem is a generalization 
of graph coloring and a special case of CSP. In that context, the central 
open problem in the complexity of the channel assignment problem is to find 
a 0*(c n )-time algorithm for a constant c independent of I or to show that 
such the algorithm does not exist, assuming ETH (or other well-established 
complexity conjecture). 
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